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Abstract. We compute the width and shape of the EPR and tunneling resonances due to dislocations in Mn- 
12 acetate crystals. Uncorrelated dislocations produce the Gaussian shape of resonances while dislocations 
bound in pairs produce the Lorentzian shape. We stress that the uniaxial spin Hamiltonian together with 
crystal defects can explain the totality of experimental data on Mn-12. 
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PACS. 75.45.+J Macroscopic quantum phenomena in magnetic systems 
nanocrystalline materials 



75.50.Tt Fine-particle systems; 



The discovery of resonant spin tunneling in Mni2 ac- 
etate |j],g] has triggered an avalanch of theoretical and 
experimental works on molecular nanomagnets (see Ref. 
H for references). Despite of a significant progress made 
in understanding Mni 2 and later discovered Feg spin-10 
systems, a number of key questions remains unanswered. 
One of them is the width and shape of the tunneling reso- 
nance. In the past, attempts were made to explain them b 
phonons [||5|] , nuclear spins g| , and dipolar fields |||7]. 
Recently, we have suggested that quantum magnetic re- 
laxation in molecular nanomagnets can be explained by 
dislocations in the crystal lattice @,|). Four recent ex- 
perimental works give evidence of the effect of defects 
on tunneling and EPR in single crystals of Mni2 and 
Fes pl| , p^ , p^|p4| , in accordance with our suggestions. The 
analysis ofthese experiments requires computation of the 
width and shape of EPR and tunneling resonances due 
to dislocations, which is done in this Letter. We argue 
that dislocations at common concentrations provide the 
observed width of tunneling resonances and the observed 
width of the EPR in Mn 12 and Fe 8 . 

Qualitatively, the importance of dislocations is clear 
from the fact that they give rise to long-ranged elastic 
strains which modulate crystal fields and thus create spa- 
tial dependence of the magnetic anisotropy. In spin tunnel- 
ing and EPR experiments the resonant values of the mag- 
netic field are determined by the anisotropy constants. For 
Mni2 crystals in the field parallel to the easy axis, with 
the Hamiltonian H = -DS 2 Z - H Z S Z + W (W being a 
small tunneling term) , the resonant spin tunneling occurs 
at 



while the EPR between the levels m and m— 1 at frequency 
lu occurs at 

H z =uj- D(2m-1). (2) 

In our two recent works |l(j,[§ we computed H' due to 
dislocations and neglected the effect of dislocations on the 
resonance condition. Meantime, the spatial dependence of 
the magnetic anisotropy D due to dislocations causes res- 
onances to spread over a certain field range. This range 
depends on the magnetoelastic coupling, the type, and 
concentration of dislocations. 

The terms in the magnetoelastic coupling that are re- 
sponsible for the modulation of the uniaxial anisotropy 
constant D can be written as 



Hu 



where 
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D' = D [g (s xx +£ y 
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1 f du a dup 
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(3) 



(4) 



is the linear deformation tensor and a, (3 — x,y,Z. The 
coupling constants go and g' must be of order one, see 
Ref. ||15| and references therein. For illustrations, we will 
use g = g'o = 1. 

For screw dislocations, one has e xx — s yy = e zz = so 
this type of dislocations does not contribute into the EPR 
and tunneling resonance linewidths. For edge dislocations 
with the axis along the z axis of the crystal and the extra 
plane z, y inserted at y > (see Fig. ||) one has u z — 
whereas other displacement components are given by [n6l 



H 7 = kD, 



0,±1,±2,...,±(2S-1) 
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Fig. 1. (a): Edge dislocation running along the y axis with the 
extra plane y, z inserted at z > 0; (b): The magnitude of the 
dislocation-induced contribution into the uniaxial anisotropy 
constant D. 



and 
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2^ 



1 -2a 

4(1 -<y) 



\n{x 2 + y 2 ) + 



2(1 -a)x 2 +y 2 



where b is the Burgers vector coinciding with the lattice 
period and < a < 1/2 is the Poisson elastic coefficient 
(we will use a = 0.25 in the numerical work). The relevant 
components of the deformation tensor are e zz = and 



(I -2a) 



II 



2tt(1 - a) x 2 + y 2 



(7) 



Displacements due to other types of edge dislocations 
can be obtained from Eqs. (g) and (J3|) by the change of 
variables. In particular, for edge dislocations with the axis 
along the y axis of the crystal and the extra plane z,y 
inserted at z > one should make a replacement y =>• z. 
This yields e yy = 0, 



and 



b (3-2a)x 2 + (l-2a)z 2 

~4ir(l-a) Z {x 2 + z 2 ) 2 



b (1 + 2a)x 2 - (1 - 2a)z 2 

47r(l-a) Z (x 2 + z 2 ) 2 



(8) 



(9) 



Generally, the axis of an edge dislocation can be di- 
rected along the x, y, and z axes of the crystal, and in 



each of these cases there are four possible orientations of 
the extra crystallographic plane. One can write 



D' = D 



9d{<P) 



(10) 



Here r is the distance from the dislocation axis, measured 
in the lattice units, whereas gD if) is a function of the 
angle which is of order one if go ~ g' ~ 1 . One can imme- 
diately see from Eq. (|l^) that the contribution of disloca- 
tions in the EPR and tunneling resonance linewidth must 
be large. Indeed, for r ~ 1 one has D' ~ D, whereas the 
spatial decay of D' is slow, so that each dislocation affects 
a large number of molecules in the crystal thus rendering 
each molecule a different value of the uniaxial anisotropy. 
This leads to a substantial inhomogeneous broadening of 
resonances which follows from Eqs. (|l|) and (0). 

In a crystal with dislocations, the deformation tensor 
at any given point is a sum of contributions due to many 
different dislocations. The superposition principle for de- 
formations follows from the linearity of the equations of 
the theory of elasticity pi and it holds everywhere out- 
side dislocation cores, i.e., for the distances from the dis- 
location axes r > 1. Statistical properties of deformations 
in a crystal depend on the spatial distribution of disloca- 
tions which is poorly known. Let us find analytically the 
distribution of the anisotropy constant D assuming that 
dislocations are distributed at random. The distribution 
function for D' in Eq. (||) can be defined as 



ffr = * 



(b'-Y^b'iv-r^ 



D' 

D' = —, (11) 
D y ' 



where ./V 3> 1 is the number of dislocations in the crystal 
and the averaging is carried out over their positions Yi 
in the plane perpendicular to the dislocation axis within 
a circular region of radius R. We choose the observation 
point in the middle of the crystal, r = 0. One can define 



N 



irR 2 irR 2 ' 



(12) 



where c is the concentration of dislocations and R c is the 
characteristic distance between dislocations. 

Let us at first analyze the large- 1 D' | asymptotes of f^, 
due to the regions with large deformations of both signs 
close to one of dislocations. In that case one can neglect 
the influence of all other dislocations and consider the one- 
dislocation model 



1 



ttR 2 Jo 



Integration yields 



dip 



rdrd D' 



fr» = 
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where D' c 
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Rc 



(13) 



(14) 



is the characteristic value of D' at the distance 
.) is the angular average. This formula becomes 
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invalid for D' <D' C , where the lines of constant D' in Eq. 
( |l3| ) cross the boundary of the region under consideration, 



r = R c . In fact, for D'<D' C Eq. (|1 

and one has to take into account otf 

( |l4| ) suggests that one should introduce the distribution 

function for the reduced quantity a 



becomes invalid 
ler dislocations. Eq. 



f a = D'J D ,, a = D'/D' c , 
which has the asymptote 

./• Q = l/H 3 , a>l. 



(15) 
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where 



irR 2 



dip / rdr exp — 



j 9d{v) 



(17) 



(18) 



In Eq. (O) we assumed for simplicity that all dislocations 
are of the same type. 

As we shall see, in Eqs. (0) and @, u ~ R c < R for 
N 3> 1, thus the argument of the exponential in Eq. ( fl8| ) is 
small and f(u>) is close to unity. Then the exponential can 
be expanded and integrated, with a log accuracy, in the 
interval \u>\ < r < R. Given that (goif)) = 0, the result 
has the form 



/M 



u 2 {g D {y) 2 ) 
R 2 



In ■ 



c R 



\uWi9Md~ 2 



(19) 



where en is a constant of order unity. Now with the use of 
Eqs. (O) and dT3) one can write 



/(w 



,A' 



w 2 (5d(<^) 2 



In- 



co-R 



exp 



~(^) 2 ln 



u\y/(9D{tp¥ 
c VN 



\w\D> 
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In the general case, with the help of the identity 2tt5(x) = 
J_ duie' luJX , the averaging over the coordinates of differ- 
ent dislocations in Eq. (0) can be factorized, 



At this point one may forget about the initial assumption 
on the circular form of the spatial region. The shape of the 
crystal only affects the value of the constant en under the 
logarithm. Eq. ( pp| ) confirms the assumption lu ~ 1/D' C ~ 
R c made above. Now we are prepared to write down the 
final result which is convenient to formulate in terms of 
the function /„ defined by Eq. ( |15| ) 



1 



du cos(cra) exp 



Vln^ . (21) 



Here the cutoff A satisfies 1 <C yl <C y/N; one cannot in- 
tegrate up to oo since the form if the integrand is only 
valid for u <§C V7V. Clearly, for large enough crystals with 
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Fig. 2. Distribution of the uniaxial anisotropy due to a ran- 
dom array of edge dislocations 



N ^> 1 the result does not depend on A. We remind that 
for the edge dislocations along the y-axis, the distribu- 
tion of transverse anisotropies is an even function. The 
distribution is shown for D' > in Fig. pi 

Integrating Eq. ( pl|) by parts three times, one can re- 
cover the asymptote of f a at \a\ ^$> 1 which is given by 
Eq. (Iq). This power-law asymptote is a consequence of 
the logarithmic singularity of the integrand in Eq. ( |2l| ) 
at u — > and it leads to the divergence of the second 
moment of f a . On the other hand, for large N the distri- 
bution function may be well approximated by Gaussian 
for not too large a. Indeed, for large N the logarithm in 
Eq. (|2l|) is weakly dependent on u and can be replaced by 
a constant. The best value of this constant corresponds to 
u for which the argument of the exponential equals one. 
This requires solving a transcendental equation that can 
be done in a perturbative way. With a good accuracy one 
can use 



In 



coVN 



=4> L = In 



c ^N\n(c VN) 



(20) which results in the aproximation 



fa 



■. exp 



a 



(22) 



(23) 



which is also shown in Fig. pi 

As the number N of dislocations in the crystal inreases, 
the function f a of Eq. ( pl| ) becomes closer and closer to 
the Gaussian, whereas the power-law asymptote given by 
Eq. ( Jig ) becomes shifted to the region of very large a 
where it is hardly visible. This effect is due to the accu- 
mulation of small contributions from dislocations situated 
at large distances from the observation point (of order of 
the linear dimension of the crystal) . Such small contribu- 
tions from distant dislocations, which lead to the Gaussian 
distribution f a , win over contributions from close disloca- 
tions responsible for Eq. (|l6| ) . Gaussian approximation for 
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the function /g,, with the help of Eq. (jla) can be written 
in the form 



fb> = 



D' 



2D 1 ,, 



: exp 



(2DL] 



where 



b'~ c = b' c VL = ^-K{g{ V T)c } 



cL. 



(24) 



(25) 



The standard deviation of D' according to Eq. ( p4| ) is 
afy, = \J2D'~. One can see that the accumulation of con- 
tributions from distant dislocations leads to the effective 
logarithmic renormalization of the concentration of dislo- 
cations c with L defined by Eq. (E3). For edge dislocations 



running along the z-axis, the quantity \/(c;_d(<£>) 2 ) is given 

by 

^MW) = -% \~ 2<T , (26) 

2V27T 1-cr 

where a is the Poisson elastic coefficient. For for g = 1 
and cr = 0.25 one has \/(gD(<p) 2 ) ~ 0.075. For edge dislo- 
cations running perpendicular to the z-axis one obtains 



V(9d(p) 2 ) = 



1 



8tt(1 - a) 
4<r(5o-ffo)( 3 5o-.9o) 



[8a 2 (50 -ffo) 2 
2 5o5o 



5 3o 2 



{9'») 2 ] l, \V) 



which for g = g' simplifies to ^{goip) 2 ) = c/ /[47r(l — 
a)]. For g = 1 and a = 0.25 one has \/(gD{p) 2 ) ~ 0.106. 

The experimentally studied Mni2 crystals are rather 
large, about 0.5 x 0.5 mm 2 , which corresponds to the cross- 
section of about 10 11 lattice cells. Even for the concentra- 
tion of dislocations as small as c = 10~ 4 per cell, the 
number of dislocation in the crystal is about N « 10 7 . 
For Co — 1 this gives L — 9.1, i.e., the effective concen- 
tration of dislocations increases by an order of magnitude, 
c = 0.91 x 10 -3 . The corresponding value of D'~ that fol- 
lows from Eqs. (|2^) and ( p7| ) for the edge dislocations 
running perpendicular to the z axis is £>- = 0.567 x 10~ 2 . 
For c = 10~ 3 one obtains L = 10.3, thus D' 5 = 1.91 x 10~ 2 . 
The renormalization of the concentration of dislocations 
and the Gaussian distribution of transverse anisotropics 
for large crystals are clearly seen in Fig. @: The distribu- 
tion broadens in the a-scale due to the increase L with N, 
Eq. H). 

The dislocation mechanism proposed in this Letter can 
qualitatively explain the experimentally observed tunnel- 
ing and EPR 00 linewidths in Mni 2 Ac. For the 
realistic concentrations of dislocations c = 10~ 3 the stan- 
dard deviation a^, = \/2D'~ rs 0.027 is in accord with the 
fit a D = Q.02D of Rcf. ||. 

Now we consider another model of distribution of dislo- 
cations in the crystal: Dislocations of opposite signs bound 
into pairs at the distance d. Distributions of this kind are 
more likely than a completely random distribution since 
here the energy of elastic strains is lower. At the distances 
r 3> d from Eq. (|l(j) one obtains 



D' 



-D^A(p), A{tp) 



d_ 
dtp 



[g D ((p) sin tp], (28) 



where ip is the angle between the vectors r and d. The 
calculations following after Eq. (11) should be now redone 
with D' given by Eq. (|28| ) and the parameters N, c, and R c 
designating the number, concentration, and the average 
distance between the dislocation pairs. For (J < iJ c < B 
the characteristic values of u> and r are u> ~ R 2 /l and 
r ~ \/uA ~ R c . The function f(u>) in Eq. (p"7J) then reads 



/(w) = 1 - 



1 



irR 2 
nl\u>\ 
~2R? 



2tt />oo 
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1 — cos 
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(29) 



This results in 



/( ^-l-*(|^)|)-exp 



Finally, Eq. 



2RI 

L7l) yields 



Trili 



2R 2 



(1^)1) 



(30) 



fb> 



1 



D' 



ir (D/)2 

irl 



(D' c f 



2RI 



(\M<p)\) = 1 ir(\M<p)\)- 



(31) 



In contrast to the random-dislocation model which is char- 
acterized by Gaussian fluctuations of the uniaxial anisotropy 
D, here the distribution of D is Lorentzian and its re- 
duced width D' c is by a factor l/R c <C 1 smaller than 
that of Eq. (|TJ). The asymptote f D , £* (1/tt)D'J(D') 2 
at D' 3> D' c is due to a single dislocation pair, as can 
be checked by an independent calculation similar to that 
for random dislocations [cf. Eq. (JL4J)]. For edge disloca- 
tions running perpendicular to the z-axis and go = g'o-, 
one has (\A(ip)\) = 3 3 / 2 5o /[47r 2 (l - cr)], which for g = 1 
and a = 0.25 yields (|A(y>)|) rs 0.175. For the concentra- 
tion of dislocations c = 10~ 3 , which is rather common, the 
average distance between dislocations is according to Eq. 
( |l2|) R c rs 17.8, in lattice units. For the size of the disloca- 
tion pair / = 5 which satisfies the applicability condition 
1<R C of Eq. (|3l]) one obtains the reduced width of the 
distribution of the anisotropy constant D' c rs 4.3 x 10 -3 
which is expectedly smaller than that for randomly dis- 
tributed single dislocations. 

If the distance between dislocations in a dislocation 
pair / is comparable with the average distance between 
dislocations R c , the distribution of D will be neither Gaus- 
sian nor Lorentzian. It still can be obtained numerically 
by the method described above. A more realistic model 
should include distribution of the dislocation-pair length 
I. We do not attempt to consider these more complicated 
models here since too little is known about the disloca- 
tions and their interaction in Mni2 and other molecular 
magnets. More detailed experimental investigation of the 
tunneling resonance and EPR lineshapes, as well as X-ray 
scattering investigations, are needed to elucidate the dis- 
tribution of dislocations in these materials. Still, the first 
results reported on in this Letter show that dislocations 
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at reasonable concentrations can be made responsible for 
the experimentally observed linewidths in Mni2. 

This work has been supported by the NSF Grant No. 
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